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SUMMARY 


Some results of studies of convergence and accuracy of finite element ap- 
proximations of certain nonlinear problems encountered in finite elasticity 
are presented. A general technique for Obtaining error bounds is also de- 
scribed together with an existence theorem. Numerical results obtained by 
solving a representative problem are also included. 


INTRODUCTION 


In this note I summarize some recent results obtained on finite element 
approximations of certain nonlinear elliptic-boundary -value problems in finite 
elasticity. The results I quote here are given in a more elaborate form else- 
where, In reference 1, Ricardo Nicolau and I reported some results on a class 
of problems in which bifurcations occur. There we consider cases in which, 
for a given set of external forces, not only can multiple solutions occur, 
but a loss of regularity can apparently result on certain solution paths. A 
complete account of these results is to be published in a lengthier article. 

The principal features of this work are (1) a priori error estimates and 
proofs of convergence of finite element approximations of highly nonlinear 
elasticity problems (these estimates are optimal), (2) errof estimates for 
multiple solutions of a nonlinear elliptic problem (these estimates are also 
optimal, but the predicated bounds are different for different solution paths), 
(3) a discussion of specific numerical results and certain special problems 
connected with the numerical analysis of this class of problems. 


NOTATION AND PRELIMINARIES 


We shall employ the following notations and conventions: 

w = (u,v,w) = displacement vector in a material body B, u, v, and w being che 
cartesian components of displacement in the material directions 
X, Y, 2. 


* This work was supported by the National Science Foundation under Grant ENG- 
75-07846. 
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Vw = gradient of w 

W - strain energy per unit volume of the body in a reference configuration, W 
being an appropriately invariant twice-continuously differentiable func- 
tion of Vw. 

V - V(w,p) = potential of the external forces per unit reference volume, p 
being a real loading parameter. 

E = 9W/9Vw = stress tensor = E(w) 

Li -.space of admissible displacements = { w: I (W + V)dXdYdZ < w = 0 on 9fi} 

(Here is a bounded open set of particles Gomposlng the interior of the 
body B and 9fi is its boundary) 

To indicate various dependences, we also use such notations as E(w), 
VV(w,p), etc. 

The potential V(w,p) is assumed to be of the form 

V(w,p) = - (pf ,w) + V^(w,p) 

where pf is a body force term and V^(w,p) is nonlinear in w. To simplify 
notations, we also introduce the operator 

f 9V 

<A(w,p),n>= I a - Vn - ~ • 11) dxdYdz ( 1 ) 

Then, formally , A is given by 

3V (w,p) 

A(w,p) = -D1 vE(w) - ' — ' (2) 

We are concerned with nonlinear bound ary- value problems of the following 
type: find w ^ U such that 


<A(w,p),n^= V n Cu (3) 

We are particularly concerned with Galerkin approximations of (3). We 
introduce a real parameter h, 0 < h ^ 1, which, of course, corresponds to the 

mesh parameter in finite element approximations, and denote {Li,} = a 

h o<h<l 

family of finite-dimensional subspaces of Li such that U is dense in U. 

0<h<l ^ 
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The Galerkin approximation of (3) then amounts to resolving the following 
problem: find w^^ £ such that 


<A(wj^,p),nh> = V ^ Bh \ 


(4) 


Upon subtracting (4) from (3) evaluated on r) = we obtain the ortho- 
gonality condition: 


<A(w,p) - A(Wj^,p),n^‘> - 0 V 


(5) 


SOME HYPGTEHSES ON THE STRESS AND POTENTIAL OPERATORS 


In many problems in finite elasticity, it appears to be justified to 
make hypotheses of the following type concerning the operator A and the space 
U: 

I. The operator A of (1) maps U into its topological dual U’ ; U is a 
reflexive Banach space with norm ] ]wj 

II. The displacement field in the body correspondin^g to a given load p 
is contained in a space U with stronger topology than (i, U being densely and 
continuously imbedded in U, 

III. The operator A is weakly continuous; i.e, if {w^} is any sequence 
converging weakly to w^, then A(w^,p) converges weakly to A(w^,p). 

IV. The operator A is coercive; i.e. 




<^A(w,p) ,w^ 

lim ’ ' ‘ ^ : = +0O 

' • 4- 1 I w l ' 


'U 


( 6 ) 


V. A sufficient condition that II holds is that A be a potential opera- 
tor with a Gateaux differential DA such that < DA(w + 0(w - w )) • ri,w - w = 0 

^ ^.o ~n -o J ~n ~o^ 

as n 00 for any sequenee {w } converging weakly to w , V n ^ U. 

VI. A sufficient condition for coerciveness is that there exists a 
constant y > 0 such that 

<A(w^,p) - ACw^.p) , w^ - ^ 2 /— ^o I 1^1 “ ^ 

where fs a positive constant and p > 1. 
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VII. There exist functions B; U x U ^ ]R and C; U x U ->1R, B weakly con- 
tinuous, such that V € U, 

Ka(wi,p) - A(w2,p),W2>1 < I 1 ^ 3 ] 1^1 - W 2 I ]^B(w^ ( 8 ) 

|<A(Wj^,p) - A(w2>p), - W2>| ^ yI ” ^2 | (9) 

where y is a positive constant and p > 0. 

Theorem 1 (Existence), Let either of the following hold; 

(i) Conditions I, III, and IV above, or 

(ii) Conditions I, IV, and V, or 

(iii) Conditions I, III, and VI, or 

(iv) Conditions I, IV, and VI. 

Then there exists at least one vector w ^ U that satisfies (3) for each 
p f € W . ■ ^ 

We emphasize that the operator A is not necessarily monotone. 


FINITE ELEMENT APPROXIMATIONS AND ERROR BOUNDS 


The subspaces ii, in (4) are assumed to be constructed using finite ele- 
n 

ment methods. Thus, the solution domain H is partitioned into E subdomains 
over which w is approximated by piecewise polynomials of degree < k. If 

^ -I - ^ 

w t U (1 U and Wj^ is its projection into Uj^ , it is well known that the subspace 
can be designed so that the following hold; 


(i) 



W, I I < C h*^j [ w 

~h"u- ° ~ 


(J 


(Id) 


h being the mesh parameter and 0 a positive number 


(ii) 




< 



V > 0 


( 11 ) 
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In (10) and (11) , and are constants independent of h. 
We proceed to determine error bounds as follows: 

1. The approximation error is e = w ” w, : 


1 1 |w 
U 

' !hl 

1 + 1 | w. - 

'U 'Uh 

■ “hi 

'u 

(by 

the 

triangle 

inequality) 

< C h'^ 
-- 0 

1 lw| 1 

+ 1 - 
U '~h 

^11 

U 

(by 

(10)) 




”h ' "hi 

'u - 


- 

IP 

i' ‘a 





(by 

(ID) 


< 

C^h^ 1 /yICA(w^ 

^,P) - 

A(w 

'Vil 

.P) . 

^h 

-Sh>l 

(by 

(9)) 


= 

c^i/Y h"! 

<'A(w,p) - 

*<"h 

;p). 

"h 

■ !h^l 

(by 

(5)) 


< 

■7 

) I Iw^ 

- w, 

1 

igi: 

W “ 

^hJ 

1 h'* 

u 

(by 

(8)) 


3. For sufficiently small h, we assume that 


B(w,w^^) ~ B(w,Wj^ - w + w) 




B(w,w) + 0 (h ) y > 0 


( 12 ) 


owing to the continuity of B (*,•). Thus 


l~h ■ ^hl h^"^^ 1 |w| I^B(w) 


(13) 


by virtue of (10), wherein B(w) = B(w',w). 

4. Combining the result 1 with (13), we see that as h 0, a positive 
constant C 2 exists such that 
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(14) 


e| I < C I |w| + h'^"'’^B(w)) 

~ ~ ^ ~ U 


Thus, for sufficiently smooth w, we obtain the optimal rate of conver- 
gence f or the nonlinear problem so long as v ^ 0 . 

Theorem 2 , Let (8), (9), and (13) hold and let there exist solutions 

to the nonlinear boundary- value problem (3). Let be a finite element 

approximation of w in a subspace in possessing properties (10) and (11). 

Then the approximation error e - w - w, satisfies the bound (14) as h 0, 

~ ~n 

Moreover, if V 0 and w is sufficiently smooth, the optimal rate of conver- 
gence is obtained for the nonlinear problem. 


AN EXAMPLE AND NUMERICAL EXPERIMENTS 

The following example is described in [1]: 

W = -E^'^nA + E^(A^ + v'^ -1) +E2(A^ + vV^ - 1)^ + E3(A^+v’^ - 2) + E^(A - 1) 

(15) 

V = - pu + K^pv^ (16) 


where A = 1 + u' (u = u(x), v = v(x)), E^,*»*,E^, are constants, and p 0, 
In this case, 

1. U-{(u,v); r (W + V)dX < <»)nwhl) 

Jg -.4 

W, (I) = Reflexive Sobolev space = { (u,v) : I (|u' | + |v' | )dx < 

Jn 


0 

u(0) = u(L) = v(0) = v(L) = O}. 

»1 


w 


u 


I o r 

wf(i) 


+ V 


kJ(I) 




*1^ dx}^ 


HI -[||u||Jl +|]v||1l ^ 

~ wjci) W^CI) 


3. a = w^(i) n wj(i) I = (o,L) 

4. p = 4j a = min(k,£-l), V = 3/2 
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The functions B(w,w) and C(w,w) are complicated functions of the compo- 
nents u and V and are given in [1]. In this case, the operator A is not mono- 
tone. 

Test problems were solved using piecewise linear finite element (k = 1). 

The problem does not have unique solutions for p > Figure 1 shows the 

computed solutions for various values of p for the case 1, 

E„ = 0.8, E„ ^ 0.5, E„ = -0.1, E/ * -0.2, K =1.0. Observe that a bifurcation 

2 3 3 4 ’ o 

is reached at p = 0.5. 

Figure 2 shows the rate of convergence actually obtained in the analysis 
computed by comparing the solution for coarse meshes with that obtained for 
100 elements. As predicted, the rate of convergence is 

0(h^ + «*= 0(h + - 0(h) 
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Figure 1.- Computed equilibrium paths. 
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